Quantum morphogenesis: A variation on Thom's catastrophe theory 
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Non-commutative propositions are characteristic of both quantum and non-quantum (sociologi- 
cal, biological, psychological) situations. In a Hilbert space model states, understood as correlations 

■ between all the possible propositions, are represented by density matrices. If systems in question 
^— ^ interact via feedback with environment their dynamics is nonlinear. Nonlinear evolutions of density 

■ matrices lead to phenomena of morphogenesis which may occur in non-commutative systems. Sev- 
' eral explicit exactly solvable models are presented, including 'birth and death of an organism' and 

. , 'development of complementary properties'. 

& : 

^ ! I. INTRODUCTION 

oo : 

René Thom's catastrophe theory is an attempt of finding a universal mathematical treatment of morphogenesis [1] 
understood as a temporally stable change of form of a system. The theory works at a meta-level and does not crucialy 
V/ depend on details of interactions that form a concrete ecosystem, organism, or society. In order to achieve this goal 
!^ the analysis must deal with qualitative classes of objects and has to possess certain universality properties. 

The purpose of the present work is similar. We define a system by an abstract space of states. The set of propositions 
y—i • which define properties of the system is in general non-Boolean. In particular, propositions corresponding to the same 
t-H , property may not be simultaneously measurable if considered at different times. Also at the same time there may 

■ exist sets of mutually inconsistent propositions. 

Although formal logical systems of this type are well known from quantum mechanics [2] it is also known that 

■ the scope of applications of non-Boolean lògic is much wider [3-7] . Practically any situation which involves contexts 
' belongs to this cathegory. Formally a context means that a logical value associated with a given proposition depcnds 

on a history of the system. In particular, an order in which qüestions are asked is not irrelevant. 
w ' The systems we shall consider are probabilistic. The morphogenesis will be described in terms of probabilities 
or uncertainties associated with given sets of propositions. The contextual nature of the propositions will require 
a representation of probabilities different from the Kolmogorovian framework [8] of sets and commuting projectors 
q-( (characteristic functions). Propositions will be represented by projectors on subspaces of a Hilbert space. 

Another element which we regard as crucial is a feedback. Feedback means that the system under consideration 
interacts with some environment. The environment is influenced by the system and the system reacts to the changes 
of the environment. Evcn simplcst models of such interactons lead effectively to nonlinear evolution equations [9]. 
Therefore, instead of modeling the interaction we will say that the feedback is present if the dynamics of the system 
^ , is nonlinear, with some restrictions on the form of nonlinearity. 

A system which interacts with environment is statistically characterized by nontrivial conditional probabilities. 
In the language of non-Kolmogorovian probability calculus this implies that states are not given by simple tensor 
produets of states. On the other hand, a simple tensor describes a state involving no correlations and hence neither 
interactions nor feedback. 

As a consequence, the nonlinearity representing feedback should disappear if the system in question and the environ- 
ment are in a product state. The latter property may be used to reduce the class of admissible nonlinear evolutions. In 
the Hilbert-space language the state of a subsystem is represented by a statistical operator p which is not a projector 
(i.e. p 2 7^ p) whenever the state of the composite system subsystem+environment is not a product state. Therefore, 
the condition p 2 — p characterizes states of subsystems which do not interact with environments. This leads to the 
following restriction: The dynamics of p is linear if p 2 = p. 

The latter condition is still not restrictive enough since it can be satisfied by both dissipative and non-dissipative 
evolutions [10] . We shall restrict the dynamics to Hamiltonian systems. In the present paper the Hamiltonian functions 
will be time independent, which roughly means that the form of the feedback does not change in time. 

Finally, we want to make the discussion universal. By this we mean two things: (1) The Hamiltonian functions 
should be typical of a very large class of dynamical systems, and (2) the results should not crucially depend on the 
form of a feedback, but more on the very fact that the feedback is present. 
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The most universal Hamiltonian functions seem to correspond to Hamiltonians with equally spaced spectra or, 
more precisely, whose spectra contain equally spaced subsets. The class of Hamiltonians includes harmònic oscillators, 
quantum fields, spin systems, ensembles of identical objects, and many others. Quite recently the rolc of Hamiltonians 
of the harmònic oscillator type was shown to be relevant to the dynamics of a stock market [11]. 

A linear Hamiltonian dynamics of p is given by the von Neumann equation 

Íp = LÒp, (1) 

with uip = [H,p]. The equation (1) may be also regarded as an abstract representation of a harmònic oscillator. 
An oscillator which occurs in many applications in biological sciences is however the nonlinear oscillator [12], whose 
abstract version reads 

ip = Çw J ·/ J ·(p). (2) 

3 

The 'genèric' equation which is the basis of our analysis is therefore the von Neumann-type equation 

ip = Ç[ffi ,&(/>)]. (3) 

3 

The index j is responsible for the possibility of having diffcrcnt parts of the system which differently interact via the 
feedback. For the sake of simplicity in this paper we restrict the analysis to only one H and a single /: 

ip=[HJ(p)]. (4) 

The only assumptions we make about / is that this is a Standard operator function in the sense accepted in spectral 
theory of self-adjoint operators, and that it should be linear whenever there is no feedback. Nontrivial examplc 
satisfying all the above requirements is an arbitrary polynomial 

f(p) = a + a 1 p+ ... + a n p n . (5) 



II. RELATION TO REACTION-DIFFUSION MODELS 



The typical reaction-diffusion models are of the form [13,14] 

iX = lüX + oj 1 f(X) (6) 

where u) = AV 2 , and A and Cj\ are, in general complex, matriecs and X, f(X) are vectors. Particular cases of (6) are 
the Swift-Hohcnbcrg, A — u, and Ginzburg-Landau models [15-18]. 

To illustrate what kind of models we arrive at consider the quadratic nonlinearity f(p) = p 2 and the harmònic 
oscillator Hamiltonian H — X)^Lo n \ n ) ( n \- ^ n ^ nc shnplest case of a one-dimensional harmònic oscillator our nonlinear 
von Neumann equation ip = [H, p 2 ] reads in position space 

ip{x, y) = (-d 2 x + d 2 y +x 2 -y 2 ) J dzp(x, z)p{z, y). (7) 

So even simplest cases lead to rather complicated integro-partial-differential nonlinear cquations. The no- feedback 
condition implies p(x,y) — ip{x)ï){y), J dz , ïp(z)tp(z) = 1, and the equation can be separated yielding the Schròdiger 
equation 

i^(x) = (-d 2 + x 2 + const)xp(x). (8) 

The Schròdinger equation may be regarded as a diffusion equation in complex time. Similarly, the nonlinear von 
Neumann equations can be mapped into diffusion type equations by replacing t by it, or by admitting non-Hcrmitian 
Hj. The Darboux techniques we are using are not restricted to Hermitian operators. 

The self-switching solutions discussed below correspond to certain p(x,y) ^ ip(x)4>{y). The 'patterns' we find in 
explícit examples are illustrated by the probability densities 

Pt,x = (x\p t \x) = p t (x,x). (9) 
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III. ENTITIES IN ENVIRONMENTS 



Considcr two Hilbert spaces: He describing an 'environment' and spanned by vectors \E), and H e describing an 
'entity' and spanned by vectors |e). The composite system 'environment+entity' is represented by either a state vector 

I*) =Y,^Ee\E,e) =Y,^Ee\E)® \é) (10) 

E,e E,e 

or by a density matrix 

P = PEE'ee>\E,e)(E',e'\ (11) 

EE'ee' 

Assuming that all expectation vàlues of random variables are represented in terms of quantum averages we can write 

(A), - (*|i4|*) (12) 

or 

(A) p =TrpA. (13) 

Of particular interest are averages representing certain statistical quantities associated only with the entities, i.e. of 
the form 

(I® A c )* = (*|J® A e |*) = Tr cPc A c (14) 

or 

(I®A e ) p = Tr p(I®A e ) = Tr c p c A c . (15) 
The reduced density matrices p c are defined, respectively, by 

Pc = Tr B |*}(*| = ]T **Ee*Ee>\e){e'\ (16) 

Eee' 

or 

p c = Tr E p = Y PEEee' \e) (e'| (17) 

Eee 1 

In particular, for product states, i.e. those of the form 

\*)=Y,^E<i>e\E,e) = W)®\4>) (18) 

E,e 

or 

P= J2 QEE«Tee>\E,e){E',e'\= g®a (19) 

EE'ee' 

the reduced density matrices are, respectively, 

p e = Tr E mm = \4>){ct>\ (20) 

and 

Pc = Tr Ep = <t. (21) 

In such a case we say that the entity is uncorrelated with the environment, i.e. probabilities of events associated with 
the entity are independent of all the events associated with the environment. 

States of composite systems are of a product form if and only if entities are uncorrelated with environments. 
Interactions of entities with environments destroy the product forms and introduce correlations. 

Reduced density matrices corresponding to nontrivial correlations satisfy the condition 
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Pc + Pe- 



(22) 



Any density matrix p is Hcrmitian and positive. From the spectral theorem it follows that there exists a basis such 
that p is diagonal. For cxamplc, any density matrix of an entity can be written in some basis as 

p c =]Tp e |e)(e| (23) 

e 

Now consider a vector 

l*)=EV^>®l e > (24) 

e 

where |\l/ e ) € TÍe are any orthonormal vectors bclonging to the Hilbert space of the environmcnt and |e) are the 
eigenvectors of p c . Then 

Tl·is|*><tt|=5>e|e><e|. (25) 

e 

In other words, for any density matrix p e one can find a state of the composite system guaranteeing that its reduced 
density matrix is identieal to p c . In what follows we shall therefore assume that a given initial p c is a result of 
correlations of the entity with the environment. If p\ p a then the corrclations are nontrivial. 



IV. FEEDBACK WITH THE ENVIRONMENT 



Typical systems discussed in the biophysics literature involve nonlincaritics given by non-polynomial functions /. 
One often encounters Hill and other functions which are continuous approximations to step functions. A simple one- 
dimensional reaction-diffusion model describing experiments on regeneration and transplantation in hydra involves 
nonlinearities with positive and negative powers [19]. The environment is here modelled by two densities describing 
concentration of activator and inhibitor producing celis. Essential to the model is the symmetry breaking of the 
two densities, a fact accounting for the nonsymmetric devclopmcnt of hydra. More refined models [20] do not need 
externally imposed inhomogeneities but involve environments acting as active chemicals. The aim of complicated 
feedback behaviors is to account for the observed symmetry breaking of the development of hydra without a need of 
putting the nonsymmetric elements by hand. 

A close quantum analogue of biophysical dynamical systems is a 'general' nonlinear von Neumann cquation (4) 
[21,22]. If / is to represent a feedback, the nonlinear effect should disappear if the entity is uncorrelated with the 
environment. Assuming the whole system is represented by a state vector the lack of correlations implies that 
\*S>) = \tp) <8> \4>) and p c = \4>){(j)\. Such a reduced density matrix satisfies pi = p c . The condition 'no correlations, no 
feedback' is formally translated into 

[H,f(p)] = [H,p] if p 2 = p. (26) 

Let us note that the above restriction means that p e = \4>){4>\ satisfies an equation which is equivalent to 

i\Í>) = H\<f>). (27) 

The latter is a general linear Schròdingcr cquation. In the absence of feedback the entity evolves according to the 
rulcs of quantum mechanics, an assumption which is rather general and weak. 

This property has also another interpretation which is entirely 'classical'. Consider a system consisting of N classical 
harmònic oscillators with freqüències uj\, . . . lun- Denote by H the diagonal matrix diag(üJi, . . . , wjv) and by \<f>) a 
column vector with entries (f>k = Çk + ipk- Then (27) is equivalent to the system of classical equations qu = túkPk, 
Pk = —üOkqk- As a consequence, the description we propose may be extended even to fully classical systems which are 
modeled by ensembles of oscillators which evolve linearly and independently in the absence of a feedback. 

Now, what are the restrictions imposed on / by (26)? As we have said before, a general density matrix has a form 
Pc = J2 e Pe\ e )( e \> wri ere p e are probabilities. The spectral theorem implies that 

f(p c ) = J2f(Pe)\e)(e\. (28) 
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The condition p\ = p c implics that p 2 e — p e whose solutions are and 1. Therefore (26) is satisfied by any / which 
fulfills /(0) = and /(l) = 1. In practical computations one can relax (26) by requiring only 

[H,f{p)]~[H,p] if P 2 =p (29) 

since having (29) one can always reparametrize the time variable t so that (26) is satisfied. The polynomial mentioned 
in the Introduction belongs to this cathegory. 

The equation (4) possesses a number of interesting general properties. For example, the quantities 

h= Tr Hf(p) (30) 
c n =Tr(p"), (31) 

for all natural n, are time independent, h is the Hamiltonian function for the dynamics and, hence, plays the role 
of the average energy of the entity (the feedback energy included). An analogous situation oceurs in nonextensive 
statistics where h has an interpretation of infernal energy [22,36]. A system with conserved h is closed. 

Conservation of c„ implies that eigenvalues of p are conserved. The latter property means that there are certain 
features of the system that oceur with time independent probabilities. However, and this is very important, the 
features themselves change in time in a way which is rather unusual in physical systems and has many analogies in 
evolution of biological systems. 



V. SOLITON MORPHOGENESIS 



There exists a class of solutions of (4) which exhibits a kind of a three-regime switching cffect [23-25] : For times 
— oo < t <C t\ the dynamics looks as if there was not feedback, then in the switching regimc t\ < t < ti a 'sudden' 
transition oceurs which drives the system into a new state which for times Í2 "C t < oo evolves again as if there 
was no feedback. Of course, the feedback is present for all times, but is 'visible' only during the switching period. 
Formally the effect is very similar to scattering between two asymptotically linear evolutions ('self-scattering'). One 
can additionally complicate the dynamics by introducing an external element which makes the form of the feedback 
time dependent. We shall illustrate the effect on explicit cxamples. 

The general equation (4) belongs to the family of cquations integrable by means of soliton methods. One begins 
with its Lax representation 

z x {i/;\ = y>\(p-\H), (32) 
-i(i,\ = ±Mf{p)- (33) 
The construction requires two additional Lax pairs 

Zv<X\ = {x\(P-vH), (34) 

-i{x\ = \{x\f{p), (35) 
z lx \tp) = {p-nH)\ip), (36) 
#) = -f{ P )W). (37) 

The method of solving (4) is based on the following theorem establishing the Darboux covariance of the Lax 
pair (32), (33) [25]. 

Theorem. Assume (V'l, (x| an d \<p) are solutions of (32), (33), (34)-(37) and Pi, are defined by 

(^IHV^I + ^P), (38) 
, 1 = (1 + ^P>(X + ^P), (39) 

p = MM. (40) 

(x\<p) 

Then 
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zx(Í>i\ = ^i\(pi-XH), (41) 

-Í{H = \{Hf(Pl): (42) 

ip\ = [H,f( Pl )}. (43) 



Let us note that the theorem is vàlid even for non-Hcrmitian H, i.c. for systems which are open. However, in the 
present paper we restrict te analysis to closed (conservative) systems chracterized by self-adjoint H. Systems whose 
average population does not change belong to this class. 

One of the strategies of finding the 'switching solutions' is the following. One begins with a seed solution p such 
that the operator 

A a := f(p) - ap, (44) 
where [a, H] = [a, p] = 0, satisfies [A a ,H] = and A a is not a múltiple of the idcntity. Now we can writc 

ip=[HJ(p)]=a[H,p] (45) 

and 

p{t) = e- iaHt p(0)e iaHt . (46) 
Taking the Lax pairs with p = v and repeating the construction from [23,24], we get 

Pl {t) = e~ ia Ht (p(0) + (ü- í/)F a (í)- 1 e- iA ^[| X (0))(x(0)|,íf]e iA ^) e iafíí , (47) 

where 

F a (í) = (x(0)|exp(i^A a í)| X (0)} 
and (x(0)| is an initial condition for the solution of the Lax pair. 

VI. 'SUDDEN' MUTATION OF POPULATION 

In our first example we consider the quadratic nonlincarity f(p) = (1 — h)p + hp 2 . The parameter h controls the 
strcngth of the fccdback. However, for any h and any density matrix satisfying p 2 = p we find f(p) = p and the 
feedback vanishes. This is consistent with our assumption that p 2 = p characterizes systems which are not interacting 
with an environment. We take the Hamiltonian H — Y^=o n \ n )( n \ which may represent a system whose energy is 
proportional to the number of its elements. Solutions of the von Neumann equation are in general infinite-dimensional 
but in order to illustrate the morphogenesis we restrict the analysis to a finite dimension. The lowest dimension where 
the effect oceurs is 3. Therefore we select a subspace spanned by three subsequent vectors \k), \k + 1), \k + 2). We 
will discuss a family, parametrized by a E R, of self-switching solutions pt = Y] m „ =0 Pmn\k + m)(k + n\ of (4). The 
solution is complctely characterized by the matrix of time-dependent coefhcients p mn . Here we only give the final 
result and postpone a detailed derivation to Sec. VIII where we analyze a generalization involving a greater number 
of 'different species'. The reader may check by a straightforward substitution that the matrix 



Poo 


Poi 


P02 


Pw 


Pll 


Pl2 


P20 


P21 


P22 



1 



15 + V5 



5 m cw \ 

m 5 + V5 £(í) (48) 

cw m 5 ) 



with 



(2 + 3i-VEi) y/3 + V5a 9e 2 ^ + (l + 4y^) a 2 

is indeed a solution of the von Neumann equation. The parameters are lüq = 1 — jg^^h, 7 = 16 ^^ h. 
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There exists a critical value ho = g 5 ^^ corresponding to lüq = 0. Using the explícit position dependence of 
the eigenstates of the harmònic oscillator Hamiltonian wc can make a plot illustrating the time dependence of the 
probability density p t , x in position space as a function of time and h. 

The dynamics we encounter in this example is particularly suggestive for h = hç, (Fig. 1) and resembles a mutation 
of the statistical ensemble described by p. The corresponding probability appears static for, roughly, — oo < t < —40 
and then also for 40 < t < oo. Switching is 'suddenly' triggered in a neighborhood of t = 0. Fig. 2 shows the evolution 
of the probability density at the origin p ti o as a function of time for diffcrent vàlues of h. For h ^ ho the probability 
density is an oscillating function of time, but in the neighborhood of t = one observes the 'mutation' which oceurs 
for any h ^ 0, the longcr the transition period the smallcr h. Duration of the switching process is of the order í/h. 
For h = the dynamics is linear (no fcedback) and there is no switching. The example shows that there oceurs a kind 
of uncertainty relation between the strength of the feedback and duration of the switching: The smaller the feedback 
the longer the switching period. 

Let us note that the probability density shown at Fig. 1 has this particular shape since we have used the position- 
space wave functions characteristic of a quantum one dimensional harmònic oscillator (a Gaussian times Hermite 
polynomials). Had we chosen any other system which is isospectral to a one-dimensional harmònic oscillator (or any 
system with equally spaced spectrum, say, a 3D harmònic oscillator) we would have obtaincd a different shape of the 
probability density. Although different choices of H imply different diffcrcntial cquations, their common feature is 
the efect of 'mutation'. 



VII. COMPOSITE ENTITIES: BIRTH AND DEATH OF AN ORGANISM 



In this example we consider an organism, that is a composite entity which undergoes the feedback process as a 
wholc. A simple model consists of a two-qubit system described by the Hamiltonian 



H = Hx <g> 1 + 1 <g> H 2 . 



(49) 



The Hamiltonian does not contain an interaction term. However, the two subentitics forming the 'organism' do not 
evolve independently. They are couplcd to each other through the feedback with the envirinment, i.e. through the 
nonlinearity. As we shall see, they become asymptotically uncoulped at t — > ±oo. In a 'distant past' the system 
consists of uncorrelated subentities which, after a period of certain joint activity, become again uncorrelated in the 
futurc. An analogy with 'birth' and 'death' is striking, and justifics the name 'organism'. 
To make the example concrete assume that 



H = 2a x C 

We will start with the non-normalized density matrix 



1 + 1 



p(o) = - 2 



which is written in such a basis that 



/5+ V7 






H = 




5- V7 
5 + 




15 



2 








1 











-1 


2 





2 


-1 



(50) 



(51) 



15/ 



(52) 



The density matrix 



p(t) = exp[-5iHt}p(0)exp[5iHt] 



(53) 



is a solution of (4) with f(p) = p 2 . Such a p(t) describes simultaneously a dynamics of two non-interacting systems 
satisfying the linear von Neumann equation 



ip = 5[2a x <E> 1 + 1 <g> a z ,p\. 
To understand why this happenes it is suficient to note that the solution satisfies 



(54) 
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[H,p 2 ] = [H,5p} = [5H,p]. 



(55) 



The environment does not trigger in this solution any switching, but only makes its evolution five times faster than 
in the absence of the feedback. The Darboux transformation when applied to p(t) produces (for more details cf. [23]) 
thc solution 



Pi(t) = exp[— 5iHt]pí nt {t) exp[5iiíí] 



(56) 



where 



Pint(í) 



5 - V7tanh2í 
_ 

13i-377-7l5+iyi05 

8cosh2t 
7t+3y7-3Vl5-i7l05 

8cosh2t 





5 + \/7tanh2í 

-15i+V7-VÏ5+i\/Ï05 
8 cosh 2J 
V7+VÏ5 
2cosh2t 



-13Í-3V7-VÏ5- 



/105 



15»- 



8cosh2t 
-Vf-VÏE-i 



8cosh2t 

/ÏStanh 2t 







-7i+3V7-3VÏ5+iVÏ05 
8 cosh 2 i 
V7+VÏ5 
2 cosh 2 i 



5-v / Ï5tanh2í ) 



(57) 



Now the switching between the two asymptotic evolutions is triggered in the neighborhood of t = 0. 

If we look at the subentities forming the organism we notice that they do not evolvc indepcndently. The casicst 
way of seeing this is to compute the reduced density matrices of the two subentities. Here we write explicitly the 
eigenvalues of the reduced density matrices. Both subsystems are two-dimensional so there are two eigenvalues for 
each reduced density matrix. They read 



p+(l) = - ± y^ÏÈ. — y^L tanh 2í, particle 1 
w 2 20 



1 ^26 + 2VÏ05 
P^=2 ± 40cosh2t ' 



The asymptotics are 



Pint(-oo) = - 



5 + V7 
5 




Pint(+00) = 



/ 5- V7 



V o 

/5 + V7 
5- V7 
5 

\ 



particle 2. 






Ï5 

5 + \/Ï5 J 

\ 


Ï5 

5- VÏ5 ) 



(58) 
(59) 



P(0), 



(60) 



(61) 



and therefore the dynamics represents asymptotically two non-interacting subentities. It is also interesting that the 
+oo asymptotics is pi(t) w p(t). At large times the 'organism' which 'dies' becomes practically indistinguishable from 
the one that never 'lived'. 

The 'life' of the organism is the period of time when the two subentities exhibit certain joint activity. Computing thc 
von Neumann entropics of reduced density matrices of the two subentities we can introduce a quantitative measure 
of this activity. Thc entropies of the two partides are shown in Fig. 3. The organism lives several units of time. 
Similar are the seales of time when thc off diagonal matrix elements of Pi nt (í) become non- negligible. It should be 
stressed that the entropy characterizing the entire organism is time independent (since eigenvalues of solutions of (4) 
are constants of motion for all /). 

Although it is clear that the 'organism' behaves during the evolution as an indivisible entity, onc should not confuse 
this indivisiblility with thc so-callcd nonseparability discussed in quantum information theory. The organism wc 
consider in the example is a two-qubit system and therefore one can check the separability of pi(t) by means of the 
Peres- Horodecki partial transposition criterion [26,27]: A two qubit density matrix p is separable if and only if its 
partial transposition is positive. It turns out that partial transposition of pi(t) is positive for any t and, hence, pi(t) 
is in this sense separable (has 'zero entanglement'). It is well known, however, that 'zero entanglement' does not 
mean 'no quantum correlations' in the system. The so called three-particle GHZ state [28] is fully cntangled at the 
three-particle level in spite of the fact that all its two-particle subsystems are described by separable density matrices. 
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VIII. MODEL WITH SEVERAL SPECIES 



The models we have considered so far corresponded to a Hilbcrt space with basis vectors \n). The only character- 
ization of a state was in terms of the quantum number n which could be regarded as the number of elements of a 
given population. Now we want to extend the description to the situation where we have a population consisting of 
several species characterized by numbers m, . . . , njv- The basis vectors are 

\n) = |m, . . . , n N ) = \m) <g> . . . <g> \n N ) (62) 

and the Hamiltonian 

H = y^(m + . . . + n N )\m, . . .,n N ){m, . . . ,n N \ (63) 

Tlj 

= J2E n \n)(n\. (64) 

n 

The Hamiltonian has equally spaced spectrum and is formally very similar to those we have encountered in the 
previous sections. The difference is that now the energy eigenstates are highly degenerated, a property which is very 
useful from the perspective of constructing multi-parameter and higher-dimensional self-switching solutions. 
For simplicity consider two species (N = 2), the quadratic nonlinearity 

f(p) = (l-h)p + hp 2 , 

and take some three energy eigenvalues Ek, Ek+ m , -E*+2m- For each energy take l + 1 vectors, which will be denoted 
by 

\0j) = \k-j,j), (65) 
\l j ) = \k + m-j,j), (66) 
|2 j ) = |fc + 2m-i,j), (67) 

j = 0, 1, . . .1 < k. We start with the unnormalized density matrix 

l 

P(0) = £>(°) (68) 

3=0 



whcrc 



Pi(0) = 1(10^X0,-1 + |2,)( 2j .|) + a+ ^ 2+ 2 4(Ò - w2) |l J )(l J l - ^%><0,| + |0 J )(2 J |). 



(69) 



Positivity of pj(0) restricts the parameters as follows: < 4m 2 < er + 46 < a . The operator 

i 

A a = p(0) 2 - ap(0) =bï- m 2 J2 \h) (h\ (70) 

3=0 



spanned by (65)-(67). We will write H = Yl\=o Hji where 



commutes with H. We denote by / and H the restrictions of the identity / and H to the 3(1 + l)-dimcnsional subspace 

- l 3- 

2 

H 3=^2(.k + nm)\n j ){n j \. (71) 

n=0 



Consider the eigenvalue problcm 



We find that the two solutions 



(^(O)-^)!^}^!^). (72) 
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V2Va 2 + 4b J ' V2 

correspond to the same j-independent eigenvalue 



vf ] ) = \h) (74) 



a+ Ja 2 + A(b-m 2 ) 

z = * ^ J - + (k + m)i (75) 

and thcrcforc 

(p(0)-iH)\<p)=z\<p). (76) 

with the same z for any 

\<p) = j2( a M 1) )+( 3 M 2) ))· ( 77 ) 

3=0 

The self-switching solution can thus be constructed by means of \ip) and reads 

Pl (t) = e-^ 1+h ^-^ Ht [p(Q) + 2iF a (t)- 1 e- hA ^)^\,H}e- hA - t y^ 1+h ^-^ Ht , (78) 

with 

F a (t) = (<p\ exp ( - 2hA a t) \v) = e~ 2hbt £ (|a,f + e 2 '"" 2 '|/3 j | 2 ) - e" 2fcW (|a| 2 + e 2 ^ 2t |/3| 2 ) . (79) 

Probabilities analogous to Fig. 1 are found if a and /i are tuned in a way which eliminates the oscillating part 
e -»(i+h(o-i))Jït ) i e for ft, = i/(i _ a ). i n this case 

p 1 (í)=p(0) + 2i(|a| 2 + e 2 - 2í /a-«)| /3 | 2 )" 1 £ [( aj |^) +^1^)) +^e^(^|), J ff] (80) 

j',j'=o 

The vectors 

!$,•(*)) - a^f) +/3,e^|^ 2) ) = |^-(í)) ® |j) (81) 

where 

i^(*)> = ^ ( - 2tm+ Sr^ ò ' m2) |fc " i} + lk+2m - J) ) + ^ lk+m - J) (82) 

are orthogonal for diffcrcnt j. Denoting H = Hi ® I + I ® H 2 one can easily compute the reduced density matrix of 
the first species, 

pf(í) = Xr 2Pl (í) = ^(0) + 2z(| Ct | 2 + e^^/C 1 — )|/3| 2 ) _1 [^ | ^- (í) ) (í ) | , ^ ] . (83) 

The entire information about the dynamics of the first species is encodcd in p[{t). Changes of properties of the species 
are given by the matrix elements 

(n\ P [(t)\n>) = (n|/ ( 0)|n') + 2i(n> - n) ^^^^ . (84) 

An immediate conclusion from the above formula is that for n — n' the expression is time independent. It follows that 
the number of elements of the ensemble does not change during the evolution. What changes are certain properties 
of the ensemble. 
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A. Example: fe = m = 1, j = 0, 1, a = 5, b = — 4 



The two-species states in the subspace in question are 



|Oo) = 


|i,o>, 


|lo} = 


|2,0), 


|2o> = 


|3,0), 


|0i) = 


|0, 1), 


|li> = 


|M>> 


|2i> = 


|2,1>- 



The two-species initial seed density matrix is given by 

5/,. „, „ „, . ,„ „, \ . 5 + V5 



Po(0) = -(|l,0)(l,0| + |3,0)(3,0|) + ^^|2,0)(2,0|--(|3,0)(l,0| + |l,0}(3,0|), 

1,1><1,1| - |(|2, 1><0, 1| + |0, 1><2, 1|), 



Pi(0) = ^(|0,l)(0,l| + |2,l)(2,l|) + ^±^ 



p(0) = p (0)+Pi(0). 
Assumc a a = ai = 1/V2, f3 = e*°/ 4 , /?i = e* 1 / 4 . Then 

1/ + 



I0o(*)> = 2(- 



1> + 



13)) 



o (*0-í)/4 



|2), 



(í)) = i(-^±^|0> + |2>) +e ^-)/ 4 |l) 



(85) 
(86) 
(87) 
(88) 
(89) 
(90) 



(91) 

(92) 
(93) 



(94) 
(95) 



Writing the restriction of H to the 6-dimcnsional subspace as 



/ 1 \ 

1 

2 

2 

3 

\0 3/ 



H 



(96) 



we can represent the two-species density matrix p\ in the form 

/ 5i 2-iVS , 



Pi = 



^1 

2 A 

2+iVE çT 
3 S 



3 , Ç 



5+vll 



3-, _j_ 2-iVE s 
2 1 - T 3 ^ 



(97) 



^1 

2 



where 1 is the 2x2 unit matrix, T denotes transposition, and 

e */4 / e ti/4 e í /4 \ 

= e t/2 + e t /2 +etl /2 ^ e ti/4 e to/4 J 

C(í) 



»*/2 



1 1 



e t/2 _|_ e to/2 _|_ e íi/2 lli 

One can verify by a straightforward calculation that (h = —j) 

i Pl = l[H, Pl }-\[H,pï}. 



(98) 
(99) 

(100) 



To illustrate the time variation of statistical quantities associated with the two-species system it is sufficient to visualizc 
the behavior of matrix elements of p\ . There are only three types of functions occuring in p\ : 
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F(t) = 
Fo(t) = 
Fi(t) = 



e t/2 


e t/2 


+ e*°/ 2 + 


etl /2 




e (í+ío)/4 




e í/2 


_|_ e ío/2 + 


e íi/2 




e (í+íi)/4 




e t/2 


+ e*°/ 2 + 


etl /2 



(101) 
(102) 
(103) 



The two parameters ío, íi, control two types of three-regime behaviors of p\. The function F is responsible for 
asymptotic properties of p\ (via (). Functions F , F\ determine properties of the switching regime (via £). For 
to <C t\ one finds Fo(t) m for all t and the switching is controlled by F(t) and í\(í); the 'moment' of switching 
is shifted proportionally to t\. For í ^> ti one finds Fi(t) ss for all t and the switching is controlled by F(t) and 
F (t); the 'moment' of switching does not depend on t\ and is determined by t - Therefore the two types of switches 
are characterized by vanishing of those matrix elements of p\ which contain either Fq or F\ . 
The asyptotic behavior of the system is given by 



F (±oo) = Fi(±oo) = F(-oo) = 0, 
F(+oo) = 1. 



(104) 
(105) 



The reduced density matrices of single species are 



p\ 



PÍ 1 



5 
2 

3 1 



2-iV5 p 
—3—^1 



2-iVÏÏ 



F 







-o" — q — + ï-^i —9 H o — ^ 



.3 + 2+iyJ^ 2+«yI Fo _ ïFi 




5+# 
-ÍF 



3 

5 
2 



2+»^ Fl _ 15+75 



(106) 



(107) 



Of coursc, all the density matrices are not normalized so that averages must be computed according to (A) 
Tr Api/ Tr pi, etc. (note that Trp! is time independent). 



IX. MORPHOGENESIS OF COMPLEMENTARITY 



According to the 'SSC theorem' [30,31] a density matrix p is uniquely determined by correlations between all the 
possible propositions associated with a given system. Each matrix element of a p can be given an interpretation in 
terms of probabilitics associated with some proposition. In the Hilbert space language a proposition is a projector, 
i.e. an operator with eigenvalues 1 and (logical 'true' and 'false'). Propositions which can be asked simultaneously 
are represented by commuting projectors. Propositions Pi, P2 which do not commute are related by an uncertainty 
relation: The more is known about P\ the less is known about P2, and vice versa. 

It is obvious that the above structures do not have to be associated with quantum systems. Just to give an example, 
many psychological tests are based on questionnaires which involve the same question asked many times in different 
contexts. The qüestions commute if the answer to a given question is always the same. However, in typical situations 
the same question has different answcrs within a single questionnaire. An ideal questionnaire involves all the possible 
qüestions asked in all the possible orders. In the Hilbert space formalism, where the qüestions are represented by 
projectors, an ideal questionnaire encodes all the possible correlations and thus, via the SSC theorem, is equivalent 
to a density matrix. 

It is also known that there exist simple examples of systems whose lògic is non-Boolean, but which do not allow 
a Hilbert space formulation [29]. The density matrix language will probably not suffice here and one has to admit a 
possibility of other state spaces and other nonlinear evolutions. The richness of available structures is immense. 

Let us finally give examples of propositions whose averages (i.e. probabilitics) change in time according to selected 
matrix elements of the self-switching solutions. The function F(t) shown in Fig. 4 is associated with the proposition 






1 


°\ 














1 











0/ 



(108) 
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corresponding to the first species as follows 

Tr Ppj(í) _ 19 + V5 + 8F(í)/3 



P(*) 



Tr/>í(í) 4 15 + v^ 



(109) 



Here p(t) is the probability of the answer 'true' associated with P. Analogously F\(t) shown at Fig. 5 is associated 
with the proposition 



Pi = 



( 1 1 
110 


Vo o o o. 



(110) 



by means of 



Pi(t) = 



TrP lP [(t) _ 1 15 + V5 + 8F!(í)/3 



Trpí(í) 



15 + V5 



(111) 



The evolution of the probabilities resembles the well known evolutions typically modcllcd by Hill functions [12] in the 
so-called sigmoidal response models [32-35,37]. Square deviations associated with the two propositions satisfy the 
uncertainty relation 



APAP 1 > - 
~ 2 



Tr [P, Pi]p{(í) VS(é' 2 + e( t+t °)/ 4 ) - (3 + >/5)e< t+tl )/ 4 



Trpí(í) 



12(15 + V5)(e*/2 + e*«/ 2 + e*i/ 2 ) 



(112) 



For any fixed t , t\, the right-hand-side of the inequality vanishes for t — > — oo and approaches v / 5/[12(15 + y/5)] for 
t — ► +oo. Fig. 6 shows this function for t\ — 0. The two propositions which were not complementary in the past 
evolve into propositions satisfying an uncertainty relation. 

In application to psychology a density matrix may represent an ideal questionnaire and, hence, a state of personality 
of a given individual. The morphogcnesis we have discussed is a simple model of development of two complementary 
concepts. The model is simplified and perhaps too far fetched. However, philosophically this is not very far from the 
approaches of Thom [1] and particularly of Zceman [38] in thcir catastrophe theory models of brain. More interesting 
in this context may be infinite dimensional cases whose preliminary analysis in terms of Darboux transformations for 
arbitrary f(p) can be found in [25]. 



X. DISCUSSION 



The model of we have described satisfies the assumptions imposed by Thom on a system of forms in evolution ( 
[1] Chapter 1.2.A). The model is continuous and the morphogenesis is a result of soliton dynamics. In this respect 
the construction is analogous to nonlinear sigmoidal response models used in biochemistry [37]. What makcs our 
construction essentially different from the models one finds in the literature is the role of non-commutativity of the 
system of propositions. 

Non-commutative propositions are related by uncertainty principies and are typical of systems which cannot, without 
an essential destruction, be separated into independent parts. The examples can be taken not only from quantum 
physics, but also from sociology (communities), psychology (personalities), or biology (organisms). In all thesc cases 
the dynamics of a system consists of two parts: One generated by infernal interactions, and the other corresponding to 
couplings with environment. We have considered only the simplcst case where the infernal dynamics is given a priori 
by a Hamiltonian of a harmònic oscillator type, and different parts of an organism (community, etc.) are coupled to 
each other only via the environment. 

The coupling with environment leads to a feedback and, hence, nonlinear evolution. The systems we consider are 
conservative, but without difficultics can be generalized to explicitly time-dependent environments or non-Hcrmition 
Hamiltonians. 

We model propositions by projectors on subspaces of a Hilbert space. States of the systems are represented by 
all the possible correlations between all the possible, even non-commuting, propositions. The choice of the Hilbert 
space language leads us therefore to a density matrix representation of states, and the dynamics is given in terms of 
nonlinear von Neumann equations. The formalism allows to consider morphogcnesis of a completely new type, for 
example a development of complementary properties. 
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The class of solutions which has an interpretation in terms of morphogenesis has fcaturcs which do not crucially 
depend on the form of the nonlinearity but more on the very presence of a feedback. The exact time development (say, 
duration) of morphogenesis does depend on initial conditions or the form of nonlinearity. However, the modifications 
do not influcnce the asymptotics, which is the qualitative element of the dynamics. For example, an organism which 
was 'born' has to 'die' but when and how will this occur depends on many details which are qualitatively irrelevant. 

Instead of conclusions let us quote René Thom's final remarks from his early work on topological models in biology 
[39]: 

"Practically any morphology can be given such a dynamical interpretation, and the choice between possible models 
may be done, frcqucntly, only by qualitative appreciacion and mathematical sense of elegance and economy. Here we 
do not deal with a scientific theory, but more precisely with a method. And this method does not lead to specific 
techniqucs, but, strictly speaking, to an art of models. What may be, in that case, the ultimate motivation to build 
such models? They satisfy, I believe, a very fundamental epistcmological need... If scientific progress is to be achieved 
by other means than pure chance and lucky guess, it rclies necessarily on a qualitative understanding of the process 
studied. Our dynamical schemes (...) provide us with a very powerful tool to reconstruct the dynamical origin of any 
morphological process. They will help us, I hope, to a better understanding of the structure of many phenomena of 
animate and inamite nature, and also, I believe, of our own structure" . 
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x 

FIG. 1. The the probability density pt, x = (x\p t \x) for the critical value /lo = (15 + 5 1//2 )/(5 + 5 1 / 2 ) as a function of time 
and x for —40 < í < 40 (in arbitrary units). The three regimes are clearly visible. The probability interpolates between 
asymptotic probabilities which are constant in time. The visible switching (morphogenesis) begins around t = —30 and takes 
approximately 30 units of time. For later times the probability density becomes indistinguishable from the new asymptotic 
state. 
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Pt,0 




FIG. 2. Probability density p t ,o at x — as a function of time and h for ho < /i < 2.45. For h > ho the switching around 
í = takes place between two different asymptotic oscillating probability densities. The switching is absent only for h = (not 
shown) where the dynamics is linear. 
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FIG. 3. Life and death of the two-qubit organism: the von Neumann entropies of partides 1 (sòlid) and 2 (dashed). The 
times where the partides are practically independent correspond to the flat parts of the plots. 
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FIG. 5. -Fi(í) as a function of t and t\ for í = 150. 
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20 

FIG. 6. Morphogenesis of complementarity. The right-hand-side of the uncertainty relation for Standard deviations AP and 
APi as a function of time t and the parameter to (íi = 0) . Propositions P and Pi are the more complementary the greater the 
value of this function. 
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